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We demonstrate that there is an optimal forcing length scale for low Prandtl number dynamo 
flows, that can significantly reduce the required energy injection rate. The investigation is based 
on simulations of the induction equation in a periodic box of size 2'kL. The flows considered are 
turbulent ABC flows forced at different forcing wavenumbers kf simulated using a subgrid turbulent 
model. The critical magnetic Reynolds number Rni^ decreases as the forcing wavenumber kf 
increases from the smallest allowed kmin = l/L. At large kf on the other hand, Rrn^ increases 
with the forcing wavenumber as Rrn^ oc in agreement with mean-field scaling prediction. At 
fc/L ~ 4 an optimal wavenumber is reached where Rmr^ obtains its minimum value. At this optimal 
wavenumber Rrn^ is smaller by more than a factor of ten than the case forced in kf — 1. This leads 
to a reduction of the energy injection rate by three orders of magnitude when compared to the case 
that the system is forced in the largest scales and thus provides a new strategy for the design of a 
fully turbulent experimental dynamo. 


Dynamo is the mechanism by which magnetic fields 
are amplified in stars and planets due to their stretch¬ 
ing by the underlying turbulent flow [ij. In the last two 
decades several experimental groups have attempted to 
reproduce the dynamo instability in the laboratory [3- 
14. The first successful dynamos were achieved in Riga 
^ and Karlsruhe 0. The flow in these dynamos were 
highly constrained and did not allow for turbulence to 
fully develop at large scales. The first fully turbulent dy¬ 
namo was achieved in where the flow was driven by 
two counter rotating propellers. However, in this experi¬ 
ment dynamo was only obtained when at least one ferro¬ 
magnetic iron propeller was used. So far other attempts 
to achieve dynamo are not successful and unconstrained 
dynamos driven just by the turbulent flows have not been 
achieved. 

One of the major difficulties to achieve liquid metal 
dynamos are the low values of magnetic Prandtl num¬ 
bers (the ratio of viscosity v to magnetic diffusivity 
77) that characterizes liquid metals and is less than 10“®. 
This implies that very large values of the Reynolds num¬ 
ber Re = UL/v (where U is the rms velocity and L is the 
domain size) are needed to reach even order one values of 
the magnetic Reynolds numbers Rm = UL/r] = Pj^Re. 
The magnetic Reynolds number is the critical parame¬ 
ter that determines the dynamo onset. For small values 
of Rm no dynamo instability exists and Rm should be 
larger than a critical value Rmc in order to generate the 
spontaneous growth of the magnetic field. The large Re 
needed to reach values of Rm above Rmc imply large 
power consumption that scales like / oc pU^L^/£f (where 
£/ is the length scale of the forcing). 


number is decreased the critical magnetic Reynolds num¬ 
ber is initially increased. The turbulent fluctuations gen¬ 
erated at large values of Re inhibit dynamo action raising 
the critical power to values much larger than in the case 
of more organized laminar flows. However when suffi¬ 
ciently large Reynolds numbers are reached this increase 
saturates and a finite value of Rmc is reached in the 
limit of Re 00 . We will refer to this value as the tur¬ 
bulent critical magnetic Reynolds number and define it 
as RmJ = lim/je-).oo Rmc- The afore mentioned stud¬ 
ies managed to reach this asymptote only by using sub¬ 
grid scale models (either hyperviscocity, a model LES, 
or a dynamical turbulent viscosity) that model the high 
Reynold number flows. It is worth pointing out that the 
different flows considered led to different values of 
implying that it is possible to optimize the flow to reduce 
Rmc ■ This is what we investigate in this work by vary¬ 
ing the length scale of the forcing if with respect to the 
domain size L. A similar study but for laminar flows was 
performed in iP. The study is based on the results of 
numerical simulations using a pseudospectral method in 


a triple periodic domain |12l Il3| and analytic estimates 


based on scale separation arguments. 

In their simplest form the governing equation for the 
evolution of the magnetic field is given by 


dth = V X (u X b) -|- rjAh 


( 1 ) 


where b is the magnetic field, and 77 the magnetic diffusiv- 
ity. u is the velocity field that is determined by solving 
the independent incompressible Navier-Stokes equation 
of a unit density p = 1 fluid. 


The dependence of the dynamo threshold Rmc on the 
Prandtl number was investigated by different groups [B“ 
for different flows with the use of numerical simula¬ 
tions. These studies showed that as the magnetic Prandtl 


dtu -I- (u • Vu) = —VP -t- i/Au -I- f. (2) 

where f is an external forcing. In the present study the 
domain considered is a triple periodic box of size 27rL 
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and f is taken to be the ABC forcing 


f = 


Asin{kfz) + C cos{kfy), 
Bsin{kfx) + Acos{kfz), 
C'sin(fc/y) + Bcos{kfx) 


with A = B = C = 1. kf is the forcing wave number, 
and we define the forcing lengthscale as if = kj^. The 
ABC flow has been the subject of many dynamo studies 
both in the laminar [ 13 , [l^ and turbulent state [ 3 ] . 

In this set up some analytical progress can be made 
in the case that kfL 1. Then one can use standard 
mean-field approximations to estimate the critical onset 


Rrn^ [ij, lla llll ■ Splitting the magnetic field in a large 


scale component B and a fluctuating part b for scale 
separation we obtain to first order for the fluctuating 
field 


r]Ah = -BVu (3) 

and for the large scale field 

^fB = V X (aB) - ?7AB (4) 

where a is a tensor such that cXijBj = (u x b)i the angu¬ 
lar brackets stand for small scale average and summation 
over the index j is implied. For the particular forcing cho¬ 
sen the a-tensor is diagonal and isotropic. The non-zero 
diagonal elements an can be calculated by solving ([3]) for 
a given velocity field and are 

— 1 U'^ 

an = —(u X A“iVjU)i = a——, (5) 

rj rjkf 

(no summation over i is implied). a is the non- 
dimensional an element that is independent of U, rj and 
kf and needs to be calculated from the flow. The growth 
rate 7 for a helical large-scale mode of wavenumber 
K = l/L is then 


ttz 

^ r -2 
7 = - VL 

■qkfL 


( 6 ) 


that leads to a critical Reynolds number Rrric = 
y^fc/L/a. This result was discussed in 0 . Note that 
this argument is true both for turbulent and laminar 
flows although the value of the coefficient a will depend 
on the level of turbulence. At large Re however a will 
reach an asymptotic value that will determine the 
value of Rm^ to be 


Rrri^ = ^ kfL/aT. (7) 

This scaling implies that in the large scale separation 
limit kfL 1 for fixed L and rj as kf is increased it 
becomes more difficult to obtain a dynamo. 

To calculate Rm^ in the absence of scale separation 
we performed numerical simulations varying the forcing 


wavenumber kf for fixed L of equations I1I2I In order 
to mimic the large Reynolds number flow that requires 
large grid size N we do not use an ordinary viscosity v 
for the dissipation but rather a dynamical wavenumber- 
dependent turbulent viscosity [l9| defined in spectral 
space as 

{k,t)= 0.27[1 + 3.58(fc/fce)*] (Ac, t)/kc (8) 

where k = |k| is the wavenumber, kc = N/3L is the max¬ 
imum wavenumber after dc-aliasing and (fc, t) is the 
kinetic energy spectrum of the flow. The same modeling 
was also used in [3| to obtain Rrn^. depends on the 
grid size N, but the dependence of the large scale com¬ 
ponents of the flow on N are expected to die off much 
faster than in the case of ordinary viscosity. 

To calculate Rm^ as a function oi kf the following 
procedure was used (see also [3,[3): For a given kf and 
grid size N a series of simulations were performed vary¬ 
ing Rm and the exponential growth rate of the magnetic 
energy was measured. The onset Rrric was determined by 
linearly interpolating the values of Rm between the slow¬ 
est growing dynamo and the slowest decaying dynamo. 
The series of runs was then repeated for higher values 
of N until either the value of Rmc remained unchanged 
in which case this determined Rm^ or we reached the 
maximum of our attainable resolution. Typically conver¬ 
gence was reached at grid sizes N = 256 but a few runs 
at N = 512 were also performed for verification. In ad¬ 
dition we solved equation [3] (with u given by equation [2]) 
with an imposed uniform magnetic field to calculate the 
elements of the a tensor and determine . 

The results for the critical Reynolds number as a func¬ 
tion of the forcing wave number kf are shown in hgure 
[1] The predicted scaling behavior © is shown with a 
dashed line in the same figure [H The black dot indicates 
the the value of calculated in Q using simulation 

of higher resolutions and a-model LES. 

The results are very motivating for future laboratory 
experiments. Although for large kfL the asymptotic scal¬ 
ing of © seems to be verified indicating that making the 
forcing length scale very small will not benefit dynamo 
experiments, at intermediate length scales Rm^ appears 
to reach a minimum around /c/L = 4 to 8. In fact the 
value of Rm^ at this optimal wavenumber is one order 
of magnitude smaller than the value of Rm^ at kfL = 1. 

The kinetic and magnetic energy spectra for slowest 
growing mode for three different forcing wave-numbers 
at highest resolution N = 512 are shown in the three 
panels of figure [H For the case that the flow is forced in 
the largest scale kfL = 1 the kinetic energy shows a clear 
^- 5/3 gpectruni while no clear power law scaling can be 
observed for the magnetic field. Most of the magnetic 
energy is concentrated in the small scales with very weak 
energy in the largest scale. 

At the other extreme where the flow is forced in the 
small scales kfL = 64 most of the kinetic energy is in the 
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FIG. 1. RtUc as a function of kfL obtained from numerical 
simulations of different resolutions N = 64—512. The data in¬ 
dicate that for N > 128, Rm^ does not vary as N is increased 
further and thus it approximates well Rrric. The error-bars 
correspond to the maximum/minimum value of Rm for which 
we obtained a clear positive/negative growth rate for the sim¬ 
ulations with N = 256. The filled circle corresponds to the 
value of RrUc obtained from direct numerical simulations in 
[^. The dashed line shows the mean field (alpha) prediction 
valid in the limit kL —>■ c». 
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small scales with a scaling in the sub-forcing scales 

and a power-law scaling for the scales larger than the 
forcing scale, suggesting equipartition of energy among 
all modes as equilibrium statistical mechanics predict 
. In addition a small peak at large scales fcL ~ 1 — 2 
is observed. The energy spectrum has been observed 
before in numerical simulations of the truncated Euler 
equations 2l| , and more recently in simulations forced in 
the small scale where the excess of energy in the largest 
scale of the system has also been observed The 

role and cause of this peak and its effect on dynamo is the 
subject of current investigations. The magnetic field on 
the other hand has a dominant peak at kL = 1, caused 
by the a dynamo that is followed by a flat spectrum 
a peak at the forcing scale and then by a power 

law until the dissipation scales. The two peaks at kL = 1 
and k = kf are in agreement with the mean field dynamo 
prediction. The two power-laws can also be explained by 
a balance between the stretching rate S of the large scale 
field by the fluctuations Ui that is proportional to 
S oc and the Ohmic dissipation that is propor¬ 
tional to Substituting oc for the turbulent 

scales and ut, oc for the scales in equipartition one 

recovers the two exponents and kP respectively. 

The spectrum was predicted in [2^ and has 

been observed in experiments 26| and numerical simula¬ 
tions [27|. The flat spectrum k° up to our knowledge is 


reported for the hrst time here. 

The case kfL = 4 that is close to the optimal wavenum- 



FIG. 2. Kinetic E{k) and magnetic Em{k) energy spectra for 
the marginally unstable modes for kfL = 1 (top), kfL = 4 
(middle), kfL = 64 (bottom). 


ber seems to be somewhere is between the two extreme 
cases. The magnetic field at the largest scale appears 
neither dominant as in the mean field case bu nor neg¬ 
ligible as in the kfL = 1 case. In the large scales there 
is not enough scale separation to observe any power-law, 
but in the small scales a power law close to k~^/^ for the 
kinetic spectrum and for the magnetic spectrum 

can be seen. 

Our results have shown that future dynamo experi¬ 
ments can benefit from forcing at scales smaller than the 
domain size by a factor of 4 to 8. To further demonstrate 
this fact in figure[3]we plot the minimum energy injection 
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FIG. 3. The non-dimensional critical injection rate In as a 
function of kjL. 

rate I = ( 27 rL)^{f • u) to achieve dynamo, normalized by 
the domain size L, the mass density p and magnetic diffu- 
sivity r] In = IL/{prf). The reason we have chosen this 
non-dimensionalization is because the mass density and 
magnetic difFusivity are properties of the liquid metals 
that vary only with temperature, while the domain size 
is typically fixed. In other words by normalizing it this 
way we ask the question: in a laboratory experiment of 
a given the domain size what is the optimal forcing scale 
to achieve dynamo with a minimal the energy injection 
rate ?. 

The result is very encouraging! The optimal injection 
rate is almost three orders of magnitude smaller than the 
case for which the forcing was in the largest scale. This 
large drop in the injection rate can be partly explained by 
considering the turbulent scaling for the energy injection 
rate I oc pU^Lr' jl. Substituting U from the definition of 
Rm we obtain / oc pRw?rf /1 and thus In oc Rm^kfL. 
Thus the energy injection rate is very sensitive on changes 
in Rm and the beneficial factor of 20 that was observed 
in figure [1] translates to a factor 2000 for the energy in¬ 
jection rate. It is also worth pointing out that the while 
for optimizing Rmf the optimal forcing wavenumber was 
between kfL = 4 and kfL = 8, in the case that In is 
optimized the optimal wave number is more clearly the 
kfL = 4. 

In the light of this result we can envision the design 
of new dynamo experiments where the flow is forced by 
an array of propellers so as to result in the small scale 
forcing required. Such experiments will challenge long 
standing theoretical assumptions of mean field dynamo 
theories. 
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